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Abstract 



D — D mixing is the source of the largest theoretical uncertainty in the extraction of 7 from B — » 
DK decays. In the Standard Model, the mixing can have a rate close to its current experimental 
upper bound and is CP-conserving to an excellent approximation. We show that neglecting CP- 
conserving D — D mixing leads to an error in the determination of 7 only at second order in the 
small parameters, Am^/^D and ATd/Vd, and is therefore very small and can be safely neglected. 
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The most precise determination of the standard model CKM phase 7 will, in the long run, 
be provided by methods based on the interference between b — > cus and b — > ucs decays 0] . 
In the case of charged B decays, the interference is between B~ — > DK~ followed by a 
£)—>•/ decay and B~ — > DK~ followed by D — > f, where / is any final state common to 
both D and D. What makes this method theoretically powerful is that there are no penguin 
contributions, and all the hadronic unknowns are in principle obtainable from experiment. 

For our purpose, it is useful to group the different methods according to the choice of 
the final state f, which can be (i) a CP- eigenstate (e.g. i^sTr , Ks4>) jJ, (ii) a flavor 
state (K + tt~) |2|, or (iii) a multi-body final state (e.g. Kstt + tt~, 7r + 7r~7r°) |3j. Additional 
variations of the basic method involve using multi-body B decays (e.g. B + — > DK + n°) 0, 
use of D°* in addition to D°, self tagging D ** states [5( and neutral B decays (both time- 
dependent and time-integrated) 0, [J. Since these measurements are statistically limited, 
an eventual combination of all the modes will be needed in order to minimize the overall 
7 measurement error 0, |^. So far, the most precise direct information on 7 comes from 
S ± — > (Ks'k + i{~)dK ± , where we use the notatio n fp to indicate a D meson decaying into 
the final state /. Both Belle jloL \v\ and BaBar jld ] have used D*° and D° decays, where 
a subtlety of a sign difference between the two D*° decay modes has been pointed out only 
recently ^J. These measurements use a sum of Breit-Wigner resonances to model the 
Dalitz plot distribution of D° — > Ks^ + ^ ■ It is possible to remove the associated modeling 
error by carrying out the measurements with a model-independent treatment of the Dalitz 

plot HQ. 

In all of the above methods, the Standard Model (SM) is assumed. (Indeed, these methods 
involve only tree-level amplitudes and, therefore, are unlikely to be affected by new physics.) 
Within the SM, the largest theoretical uncertainty is due to D — D mixing. The parameters 
that describe the mixing are 

x _AmD^ y = ^£R^ = arg(-), (1) 

where Ar^ (Amo) is the decay width (mass) difference between the two neutral D mass 
eigenstates, Y D {m D ) is the average decay width (mass) of the mass eigenstates, and q 
and p are the elements of the rotation matrix between the interaction and mass eigenstate 
bases 0. (Here, and in what follows, we choose a phase convention such that the tree-level 
D decay amplitudes are real, see Eq. (J2J)). The parameters x and y cannot be calculated 
reliably in the SM. In particular, one cannot rule out the possibility that they are as large 
as x ~ y ~ 0(1CT 2 ) [la], which is the range experiments are beginning to probe 13, Il8| . 
A robust SM prediction, however, is that D — D mixing is CP-conserving to a very high 
accuracy, with a CP- violating phase of order 6 ~ O(10 -4 ) jl8| . 

If the D — D mixing parameters are known, their effect can be corrected for in the 
measurement of 7 0, Eo|- Without knowing their values, assuming x — y — introduces 
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an error in the extracted value of 7. Making that assumption can introduce an error in 
the determination of the branching ratios used in the ADS method [2] of the order of 
x/rf ~ y/rf ^ 20%, where 77 is defined in Eq. (J2J) below. Naively, one may conclude 
that this introduces a similar error in the extracted value of 7. It is the purpose of this note 
to explain why this is not the case. We find that the effect is at most quadratic in x and y, 
thus it is very small and can be safely neglected. 

Let us review the approach of extracting 7 neglecting D — D mixing. We choose the 
phase convention in which the D meson decay amplitudes have the form 

A(D° -> f) = A f , A(D° -+f)=A, = A f r f e- iS f, (2) 

such that Af and 77 are positive. Since in the SM there is essentially no CP violation in the 
D decays, 5f is a strong phase difference. Due to the abundance of flavor-tagged D decays 
at the B factories, the values of Af and 77 can measured very precisely from the decay rates 

r / = T(D° - /) = A% Tf = T(D° -> /) = A)r). (3) 
In the absence of D — D mixing, the amplitude for the cascade decay B + — > f B K + is 

A(B + - f D K + ) = A B A f [r f e- l5 f + r B e^ B+ ^] , (4) 

where A B and r B are positive parameters, 5 B is a strong phase difference, and we have 
defined 

A B = A{B + -> D°K + ), A B r B e l{&B+l) = A{B + -> D°K + ). (5) 
The sensitivity to 7 comes from the interference term in the decay width 

T(B + -> f D K + ) = A\A)\r) + r| + 2r B r f cos(5 B + 7 + 5 f )\. (6) 

A similar expression for the B~ — > JdK~ decay width is obtained by making use of the 
absence of direct CP violation in the D decay: 

T(B- -> f D K-) = A\A)\r) + r| + 2r B r f cos(5 B - 7 + S f )). (7) 

Each final state / introduces two new observables, T(B~ — > fr>K~) and — > f B K + ), 
and a single new unknown, 5/. The unknowns describing the i? — > DK part of the cascade 
decay, r B , 5 B , and 7, are the same for all D decay final states. Therefore, with enough D 
decay modes, there are more observables than unknowns, and the values of all the unknowns 
can be determined. 

We now study the effect of CP-conserving D — D mixing. Specifically, we ask what error, 
A7, is introduced in the extracted value of 7 when the analysis is done assuming no D — D 
mixing. A crucial ingredient in the approach described above for extracting 7 is that 5f 
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is a pure phase, i.e., a single real parameter. Specifically, we assumed that the absolute 
magnitude of the interference term 



\A f A}\ = \Ajr f e^\ = Ajr f , (8) 

which is used in Eqs.dHJ) and (J7J), is already measured in flavor-tagged D decays. However, 
due to D — D mixing, time evolution dilutes the absolute magnitude of the interference 
term (JHJ), which becomes another unknown. It is the deviation of the magnitude of the 
interference term from its naive value that introduces the error in the extracted value of 7. 



In t 
are [l5| 



re presence of CP-conserving D — D mixing, the time-dependent D decay amplitudes 

A f (t) = A(D°(t) - /) = g + (t)A f + g_(t)A fi 
A f (t) = A(D°(t) - /) = g + (t)A f + g4t)A ft 
where the time evolution functions are 

g+{t) = exp(— imz)t — r/2) [cosh (yr/2) cos (xr/2) + i sinh (yr/2) sin (xr/2)] 

w exp(-im D t - r/2) [l + (y + ixf t 2 /A\ , 
g~(t) = exp(— imz)t — r/2) [— sinh (yr/2) cos (xt/2) — i cosh (yr/2) sin (xr/2)] 

w exp(—im D t — r/2) [(— ia; — y) r/2] , (10) 

with r = r^t. The approximations in Eqs. ()10j) hold to second order in x and y. D — D 
mixing changes the time- integrated decay rates of Eq. (jSJ) at leading order in x and y. The 
precise change is not important for our purpose. The key point is that the time-integrated 
decay rates 

T f = [ dt\A f (t)\ 2 , T f = I dt\A f (t)\ 2 , (11) 



which are measured in tagged D decays, are exactly the rates that enter the B decay rate, 
which now reads 



T{B + - f D K + ) = A 



12) 



f / + r 2 B T f + 2r B Re (V (<5s+7) J dtAf(t)A f (t)* 
The impact of D — D mixing on the 7 measurement occurs only in the interference term 

J dtA f (t)A f (ty = y/rfF f e a 'e*', (13) 
where 5f is a pure strong phase 1 and 

e f = ~(x 2 + y 2 ) (js+nj -\{x 2 cos26 f + y 2 sin26 f ) (14) 



1 Note that this is not the case in the presence of CP-violating D — D mixing when the phase Sf is a 
combination of a weak and a strong phase. 
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describes the dilution due to D — D mixing. The parameter e/ gives the approximate mag- 
nitude of the shift A7 in the determination of 7. Since the leading term in e/ is proportional 
to (x 2 + y 2 )/r 2 , A7 is larger for cases where rf is smaller. Apart from the trivial case of no 
mixing (x = y = 0), e/ vanishes only if r/ = 1 and either y = and 8 = kir, or x = and 
5 = 7r/2 + fc7r, where k is an integer. In all other cases, 6f is positive. 

In the special case e/ = 0, there is no change in the 7 measurement. Each new mode / still 
introduces only one new parameter, 8f, which is obtained from the fit to the decay widths. 
Moreover, the form of the equations is unchanged, as can be seen by defining Tj = Aj, 
Yf — Ajr 2 and comparing (fTTj) . (fT2*j) with ©• Therefore, for 6f — 0, the correct value of 7 
is measured, even if the formalism used in the analysis ignores D — D mixing. 

Our first main point is that while in general tf 7^ 0, tf is of second order in the small 
parameters x and y. Therefore, its effect on the measurement of 7 is small. Moreover, given 
measurements of, or upper limits on, x and y, the impact of D — D mixing can be accounted 
for without the need to perform a time-dependent analysis of the B decay. This can be done 
by using (ITU) and (IT31) in (fT2l. 

That €f is of second order in x and y can be understood as follows. One can think of 
integration over time as a scalar product in the vector space of time-dependent complex 
functions 

(A f ,A,) = J dtA f (t)A f (ty . (15) 

Then 

\(A f ,Af)\ 2 = (A f ,A f )(Af,Af)\ cos A| 2 , (16) 
where A is a small angle linear in x, y. The difference that defines ej in Eq. (|TH|) is then 

e f oc (A f ,A f ){Af,Af) - \(A f ,A f )\ 2 oc 1 - cos 2 A ~ 0(A 2 ) ~ 0(x 2 ,y 2 ). (17) 

We provide an explicit expression for A7 in one specific case where 7 is extracted from 
a combination of a doubly Cabibbo-suppressed decay width T(i? ± — > (K T 7i :iz )r)K ± ) and a 
decay width r(5 ± — > {fcp)nK ± ) into a CP eigenstate. To first order in and 77 (here 
/ = Kir), we get 

^ cos 7 sin 27 

^ cos7[cos2(<5/ + 5b) — cos25 B ] + {r B /rj) cos(<5/ + <5b)[cos 27 — cos 25b] 

where e/ is taken from (^H). For x 2 + y 2 ~ 2%, 7 ~ 60°, re ~ 0.2, and rj ~ 6% we find the 
typical range A7 ~ 0.1 — 1°, depending on the values of the strong phases Sf and 8b- By the 
time the precision of the 7 measurement reaches this level, we will have either measurements 
or tighter upper limits on x and y, so that the measurement could be corrected for this shift. 

Next, we consider the effect of D — D mixing in the case of multi-body D decays. For the 
Breit-Wigner treatment of Dalitz plot, the corrections due to D — D mixing arise at 0(x 2 , y 2 ) 
as in the two-body case discussed above. Similar considerations apply in both the two-body 
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and three-body cases, with the difference being that in the Breit-Wigner-based Dalitz plot 
analysis, 77 varies over the Dalitz plot. The lowest value of 77, and hence the largest A7, 
is obtained in areas populated by doubly Cabibbo-suppressed decays. Specifically, for the 
final state / = Kstt + tt~ this is the region of the decay D° — > K* + n~ , which contributes 
most to the 7 measurement. Nonetheless, with 77 of order a few percent in this region, 
this still results in a small contribution to A7. Moreover, the overall value of A7 is smaller 
due to the contributions of other regions in the Dalitz plot in which 77 is larger. The 
shift A7 is significantly smaller for singly Cabibbo-suppressed multi-body decays, in which 
77 ~ 0(1) 0- 

Our second main point is that CP-conserving D — D mixing does not affect the determi- 
nation of 7 if the relevant Dalitz plot parameters are determined by binning the Dalitz plot 
according to the model-independent approach of Ref jj. The phase space integration over 
bin i of the Dalitz plot introduces two new real variables, q and sf 

Ci + isi = ° l lSl , (19) 

J- i 

where 

Ci + isi = J dp J dtAf(t)A f (t)*, Tt = J dp J dt\A f {t)\ 2 . (20) 

The variables q and s« are determined either from the binned Dalitz plot obtained from the 
B decay sample, or from time-integrated decays of entangled D states at a charm factory 
operating at the \&(3770) [3J. The point is that measuring q and s, already accounts for 
the dilution due to D — D mixing. This is demonstrated by the fact that in the two-body 
case, one can replace the two variables 5f and e/ of (fT3|) with q and s,, which satisfy 
c 2 + s 2 = 1 — 0(x 2 ,y 2 ). The method of Ref. 0| is already designed to handle c 2 + s 2 < 1, 
which in multi-body decays arises due to the phase space integration over each bin. 

We concentrated on the case of CP-conserving D — D mixing, since this is the case in 
the SM. With new physics, this may not be the case. Then, our results do not hold and 
larger effects are introduced. For example, consider the case where there is new physics in 
the mixing, with a CP- violating phase 9 ~ 0(1). Then the assumption of no D — D mixing 
introduces an error in the value of 7 of order A7 ~ 0(x6,y6) which is linear in the small 
parameters. 

It is instructive to compare our results to those of 0. The analysis in Q corresponds 
to a situation in which the D decay amplitudes were determined from D decay data by 
taking D — D mixing into account, but neglecting mixing in the B decay analysis. The error 
in the value of 7 extracted in this way is linear in x and y, regardless of whether D — D 
mixing is CP-conserving or not. Here, on the other hand, we show that when both the D 
and the B decay amplitudes are extracted ignoring D — D mixing, CP conserving D — D 
mixing induces an error in the extracted value of 7 that is of second order in x and y. This 
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provides a simpler practical approach for solving the problem introduced by CP-conserving 
D — D mixing, which is the case in the SM. 

To conclude, we show that within the SM, neglecting D — D mixing in the extraction of 
7 using B — > DK type decays introduces at most an 0(x 2 ,y 2 ) effect. This is a very small 
effect that can be neglected for all practical purposes. 
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